Bound-free pair production cross section in heavy-ion colliders 
from the equivalent photon approach 
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Abstract. - Exact calculations of the electron-positron pair production by a single photon in the 
Coulomb field of a nucleus with simultaneous capture of the electron into the K-shell are discussed 
for different nuclear charges. Using the equivalent photon method of Weizsacker and Williams, a 
simple expression for the bound-free production of e^e~ pairs by colliding very- high-energy fully 
stripped heavy ions is derived for nuclei of arbitrary charge. 
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Introduction. — There has been some theoretical and 
experimental interest in the process of single-photon pair 
production during the last few decades. Precise knowledge 
of the single-photon bound-free pair production cross sec- 
tion in the electrostatic field of a fully stripped nucleus 
allows to calculate the bound-free pair production at the 
ipteraction point of a heavy ion collider by using the equiv- 
alent photon method of Weizsacker and Williams [1-3]. 
lilecently, first observations of beam losses have been re- 
ported from measurements performed at the BNL Rel- 
ativistic Heavy Ion Collider (RHIC), where a beam of 
^^Cu^^"'' ions with 100 GeV/nucleon has been used [4]. 
Xn this paper, we derive a simple and compact expression 
for the bound-free pair production in relativistic heavy- 
ion collisions, which is expected to be a major luminosity 
limit for the Large Hadron Collider (LHC) when it oper- 
ates with heavy ions because the localized energy deposi- 
tion by the lost ions which have captured an electron may 
quench superconducting magnet coils [5]. 

The equivalent photon method. The simplest 
way to get a reliable estimate of the cross section in RHI 
collisions is provided by the equivalent photon method, 
which is originally due to Fermi [6], and later on devel- 
oped by Weizsacker and Williams. Within this frame- 
work, the target ion is considered as fixed. The projec- 
tile with impact parameter b is assumed to move on a 
straight line with velocity v and a relativistic Lorentz fac- 
tor 7p = 27^ — 1, where 7c is the Lorentz boost of the ions 
in the center-of-mass frame. The projectile is accompanied 
by its contracted electromagnetic field, which corresponds 



to a spectrum of equivalent photons, given by 
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where Kq{Ki) arc the modified Bessel functions of the 
second kind of order zero (one), x ~ ujb/^pV, and a is the 
fine-structure constant. The cross section for an electro- 
magnetic process in a highly relativistic collision is then 
obtained by integrating the single-quantum cross section 
over the frequency spectrum and from a minimum impact 
parameter 6„im = R^ which was, in our case, chosen to 
be the Compton wavelength as the typical length scale for 
pair production, to infinity: 
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where a is now the total cross section of the electromag- 
netic process. The integration of N{lo, b) over b can be 
carried out to give [7] 
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where C = uiR/^pV is an adiabatic cutoff parameter. For 
7^1 (except for extreme low-energy frequencies, satisfy- 
ing the relationship ljR/c <C 1), one can use the excellent 
approximation 



X 10 



n{ui) = -Zlalog 

TT 



-^,^alog(^- 



(5) 



where 5 = 0.68108... is a number related to Euler's 
constant. In the limit of very large frequencies, oj 3> 
jpv/R, an adiabatic cutoff sets in and one has n(w) = 
(aZp/2)e~^'^^/'''p^. For practical calculations, it is usu- 
ally sufficient to set the velocity v equal to the speed of 
light c. 

Calculations of the single-photon bound-free pair pro- 
duction cross section cr^g_ with simultaneous capture of 
the electron into the K-shell of a target nucleus with nu- 
clear charge number Zt have been presented in [8]. In 
the meantime, the program used in [8] has been refined to 
higher precision, leading to small corrections of the orig- 
inally calculated K-shell cross sections of some few per- 
cent. In [9], it has been shown that the higher shells 
contribute about 20% of the K-shell capture to the to- 
tal bound-free pair production cross section cTg+g-, i.e., 
CTe+e- — C(3)o'^g- provides a good approximation for 
(7e+e-, where is the Riemann zeta function defined by 
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Fig. 1: Single-photon cross section for pair production with 
K-shell capture as a function of tlie photon energy e-y — Tiuj 
for different nuclear charges Zt. The cross sections are scaled 
with Zf. The dash-dotted curve displays the Sauter cross sec- 
tion according to eq. (8), which agrees well with the numerical 
result for Zt = 1 (upper solid curve). The dashed curve corre- 
sponds to Zt — 29 (copper), the lower solid curve to Zt = 92 
(uranium). The dotted curve shows the high-energy limit of 
the Sauter cross section according to eq. (7). 
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for I > 1, and C(3) ~ 1.2020569. 



Single-photon cross sections. — Theoretical cal- 
culations using semirelativistic Sommerfeld-Maue wave 
functions for the positron and a semi-relativistic approx- 
imation for the bound electron wave functions valid for 
aZt <C 1 show that, for Tico ^ 1 and light nuclei, the cross 
section ct^ _ is approximately given by the high-energy 
limit of the so-called Sauter cross section [10] 



(7) 



where me is the electron mass and Ac the Compton wave- 
length. 

The high-energy limit of the Sauter cross section over- 
estimates the exact cross sections at low photon energies. 
Indeed, Sauter originally derived an exact result in the 
limit of vanishing Zt for the K-shell photoelectric effect, 
which can be transcribed into the corresponding cross sec- 
tion for pair creation [11] (me ^ h ~ c = I) 
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where and p+ = (E^ — 1)^/^ is the energy and the mo- 
mentum of the positron, respectively. The photon thresh- 
old energy for pair production ujth = 1 + 7e is given 
by the sum of the positron rest energy and the energy 
7e = (1 — a'^Zfy^'^ of the bound electron in the K-shcU. 
Generally, the positron energy is given by E^ = lu — "fe, 
and in the limit Zt one obviously has to = £'+ -I- 1. 

A comparison of the Sauter cross section and exact nu- 
merical calculations for Zt = 1, 29, and 92, is shown in 
Fig. 1. The numerical result for Zt = 1 agrees well with 
eq. (8). However, at high energies, the cross section for 
hydrogen is smaller by a factor of 0.971 than the Sauter 
cross section. This effect becomes more pronounced for 
target nuclei with large Zt, such that eq. (7) must be 
modified by an additional factor f{Zt) according to 
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Some values for [Zt, f{Zt)] obtained from exact numeri- 
cal calculations are [92,0.196], [82,0.216], [29,0.484] and 
[1,0.971]. These values agree very well with those found 
by Pratt [12] and coincide with the values calculated by 
Agger and Sorcnscn [13]. 

A useful approximation of f{Zt) is given by 
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j-g-j which reproduces f{Zt) to an accuracy of 2% in the phys- 
ically relevant range < Zt < 92. 
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As a further effect one observes that the position and the 
shape of the cross section peak varies for different nuclear 
charges. Whereas the peak is shifted to higher photon 
energies for hght and medium-heavy nuclei, it is located 
at lower energies with respect to the Sauter cross section 
peak for > 62. 

Pair production in heavy-ion collisions. Eq. 

(8) in conjunction with the equivalent photon method al- 
lows to derive a precise numerical expression for the total 
bound-free pair production in relativistic heavy-ion colli- 
sions in the case of small target charge aZt ^ 1 (with 
R = 1 and w = 1 in eq. (5)), given by 
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with ci = 0.21746 and C2 = 0.58706. Eq. (11) stiU over- 
estimates the cross section, since the single-photon cross 
section is suppressed by a factor f{Zt) at large photon en- 
ergies. Additionally, the shape of the single photon cross 
section peak has some impact on af/ni, since the photon 
flux generated by the Coulomb field of the projectile be- 
comes large at low equivalent photon energies. The defect 
can be remedied by a minor modification of eq. (11) 
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where a numerical fit for < Zf < 92 gives 



The coefficients A and B were obtained from numerical 
calculations. In our case, the coefficients follow directly 
from eqns. (13,14,15). 

Cross sections for several nuclei and energies have been 
calculated, e.g., by Meier et al. [9] and Baltz et al. [15]. 
For a RHIC Cu-Cu collision with 100 GeV/nucleon, one 
obtains with 7c 107 and 7^ = 27^ — 1 a cross section of 
0.198 barn, in agreement with recent experimental mea- 
surements [4]. In the case of the LHC, one obtains from eq. 
(15) for a Pb-Pb collision with 7c = 2957 a cross section of 
220 barn for K-shell capture, compared to 225 barn in [9]. 
For Zp = Zt = 20 and 7c = 125 (7c = 3750), the cross 
section for K-shell capture obtained from the equivalent 
photon method is 1.57- 10"^ b^j.^ (2.96- 10"^ barn). The 
corresponding values found in [9] are 1.61 ■ 10~^ barn and 
2.92 • 10~^ barn, respectively. A cross section of 89 barn 
for Au-Au collisions at RHIC energy given in [15] is also 
in excellent agreement with eq. (15). This demonstrates 
the validity of the equivalent photon method for the full 
physical range of charge numbers Z and large 7c. 

As an interesting final remark, we briefly compare the 
results obtained above to approaches presented in [7] 
and in [16]. Based on Sommerfeld-Maue wave functions 
for positrons and semi-relativistic electron wave functions 
constructed from the non-relativistic electron wave func- 
tions by including perturbative correction terms of the or- 
der aZ^ an approximate expression for pair production 
with K-shell capture was derived in [7] {Z = Zt = Zp) 



33^2 7-77 27raZ 
— \,.a Z - 
20 ^ 



1 



\ogiSjp/2) - 5/3 . (17) 



The numerical factor 33/20 — 1.65 agrees in a satisfactory 
manner with the corresponding factor 8ci = 1.74 appear- 
ing in eq. (11). The deviation of the cross section from 
the Z^-scaling is basically described by 
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with ci^2 = ci,2(0) given above. 

Expressing the cross sections directly in barn leads to 

aBHi = 1.58 ■ 10-^' Z'^Z^f{Zt)[ci{Zt) log(7p) - C2{Zt)]. 

(15) 

The formula above describes the equivalent photon cross 
section with an accuracy better than 2% for Zt < 92 and 
typical RHIC and LHC energies. It should be mentioned 
that it has already been suggested in [14] that the cross 
section anni can be established in the form 



ctrhi = ^log7p + B. 



(16) 



A comparison of f{Z) with f{Z) obtained from exact nu- 
merical calculations is presented in Fig. 2. Whereas f{Z) 
is very acceptable for aZ ^ 1, it leads to an underestima- 
tion of cross sections for heavy nuclei by a factor of about 
2. The approach presented in [16] uses the same bound 
electron wave functions as in [7], however, positron wave 
functions are described by normalized plane waves in con- 
junction with a correction term proportional to aZ . The 
approximation in [16] suffers from the opposite effect, i.e., 
cross sections for nuclei with large charge numbers turn 
out to be too large by factors exceeding 2. This observa- 
tion shows that is it crucial to work with exact solutions 
of the Dirac equation in order to get reliable results in the 
present case. 

Conclusions. — A heuristic analytic fit for the char- 
acteristic function f{Z) describing the single-photon pair 
production with capture on fully stripped nuclei at high 
energies is presented in eq. (10). This fit allows to describe 
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Fig. 2: Comparison of the Z-dependence of o^^^^, expressed 
by f{Z) from eq. (18) (solid line), and exact values for f[Z) 
(dashed line). 



the single-photon pair production cross section at high en- 
ergies with an accuracy comparable to presently available 
numerical calculations. Furthermore, eq. (15) in conjunc- 
tion with eqns. (10,13,14) allows to evaluate the bound- 
free pair production cross section in relativistic heavy-ion 
collisions and the corresponding coefficients A and B in 
eq. (16) in a simple and reliable manner. It should be 
mentioned that in currently available calculations, multi- 
photon interactions between the photon-emitting nuclei 
are neglected, which might lead to an additional correc- 
tion of the cross section of a few percent, probably with 
different energy dependence from the overall cross-section. 
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